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ABSTRACT 



We analyse the properties of the host galaxy of the optically dark gamma-ray burst (GRB) 020819 (z = 0.41) and discuss the possible 
implications in the context of "dark" GRBs. We present g'r'i'z'JHK photometry of the host galaxy and fit the broad spectral energy 
distribution including the public Spitzer IRAC data using stellar population models. The broad spectral energy distribution (SED) 
indicates a high extinction. Ay ~1.8 - 2.6 mag, for this relatively massive galaxy. This is the highest global extinction for a GRB host 
galaxy with a robust spectroscopic redshift. The properties of the host galaxy are indicative of dusty, intense star-formation, which 
differ from those of the current sample of GRB hosts. This implies that the dust extinction is one of the main reasons for the darkness 
of low-redshift bursts and that the long GRB host population is far more diverse than previously anticipated. 
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1. Introduction 

The afterglows of gamma-ray bursts (GRBs) were predicted 
by the canonical fireball theory (e.g.. iMeszaros & Rees |[l997h 
and have since been detected for more than a decade 
(Ivan Paradiis et al.1 Il997h . The optical-to-X-ray flux ratio of 
a GRB afterglo w is determined by the fireball theory (e.g., 
ISariet al7lll998l) . where the optical-to-X-ray spectral index of 
a canonical GRB afterglow is expected to be fiox ^ 0.5 
dJakobsson et al. 2004). However, for a significant proportion 
(25% - 50%) of all well-localised GRBs, no optical/near-infrared 
afterglow is detected and/or the optical afterglow emission 
is lower than that exp ected from the X- r ay afterglow emis- 
sion (pnv < 0.5) (e.g.. iFvnboetal. 1 120091: ICenkoetani2009t 
Ide Pasquale et al. I l2006h. These bursts are called "dark" GRBs. 
The nature of the dark bursts is still to be understood, although 
several ideas have been proposed to explain why some bursts are 
dark in the optical bands: (i) shifting of the rest-frame optical 
afterglow emission and the Lyman-limit towards longer wave- 
lengths for high redshift bursts, (ii) intrinsic dimness of the after- 
glow (for those with no optical emission and unknown Pox), or 
(iii) high extinction in either the host galaxy or the circumburst 
environment and along the line of sight (e.g. , Fynbo et al. 2001; 
iLazzati et al. ll2002Llde Pasquale et al. ll2003LlFox et al. Il2003i) . 

GRB 0208 ISQ was detected by the High Energy Transient 
Explorer (HETE) satellite and identified as a moderately bright 
long-soft burst with T 90 ~20 seconds and a peak brightness 
of ~5 crab ( Vanderspek et al. 1 12002) . Follow-up ground-based 
observations led to a /T'-band limit of 19.5 mag and R-band 
limit of 22.15 for the afterglow emission ~9 hrs after the burst 



' This burst is in fact GRB 0208 19B (see Jochen Greiner's web site 
Ihttp://www.mpe.mpg.de/~jcg/grb020819.htmll. However it is referred 
to as GRB 020819 in all related GCNs and published papers, and we 
therefore also refer to it as GRB 020819 to avoid further confusion. 



(iKlose et al. ll2"003l:lLevan et al. ll2002h . Therefore, GRB 020819 
wa s classified as a dark bu rst. 

iFrail & Berger I (l2002h later discovered a declining radio 
source consistent with radio afterglow emission located at a faint 
blob 3" away from the centre of a spiral galaxy CJakobsson et al. I 
2005). This galaxy, at z=0.41, is identified as the host galaxy 
of GRB 020819 with a chance coincidence probability of 0.8% 
( Jakobsson et al. 2005). The host has a visible extent of ~6"- 
7"(32.5 - 38 kpc) in the /?-band (see Fig.[B. 

Here we present broad-band photometric observations of the 
host galaxy of GRB 020819 (Sec. 2) and the analysis of its spec- 
tral energy distribution (Sect. 3). In Sect. 4, we discuss the results 
and their implications for the dark GRB framework. 



2. Observations 

The host galaxy of GRB 020819 was observed with GROND 
(lGreineretal.112008 ') in g'r'i'z'JHK on 29 October 2007 for ~1 
hour at an average airmass of 1.23. The data reduction and aper- 
ture photometry of the host was performed using the data re- 
duction and p hotometry to ols of the GROND pipeline (see e.g. 
Kupcu Yoldas et al.1l2008h based on IRAF/PyRAF Astrometry 
and photometric c alibration wa s performed using the SDSS cat- 
alogue ( Abazaiian et al. Il2 009[) for the o ptical bands, and the 
2MASS catalogue ( Ski-utski e"et al. Il2006h for infrared bands. 

The brightness of the host galaxy was derived using aper- 
ture photometry with a radius of ~4.5"(3xPSF) in all bands. 
The results of the photometry are shown in Table [1] All mag- 
nitudes are corrected f or a foreg round Galactic extinction of 
E(B-V) = 0.069 ( Schle geletal. 119 98). We note that our g'-band 
magnitude, 20.31+0.02, is brighter than the B-band magnitude, 
21.9+0.5, of the host reported by Jakobsson et al. (2005). The 
difference is larger than that expected because of the differences 
between the two filter curves. The reason for this difference 
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Table 1. The photometry of the host galaxy of GRB 0208 19 " 



J 



H 



K 



20.31±0.02 19.43±0.06 19.10±0.07 18.73±0.01 18.46±0.04 18.23±0.04 18.46±0.11 



" All magnitudes are in units of AB mag and are corrected for foreground Galactic extinction. 



is pro bably the inaccur a te ca libration of the Keck B-band data 
where 'Jakobsson et al. I (1200 5) used a single faint star for cali- 
bration, which was the only available star in the field (Jakobsson, 
private communication). 




Fig. 1. VLT/FORS2 /?-band image of the host galaxy of GRB 
020819 taken on 15 Sep 2002. It is clearly resolved as a barred 
spiral galaxy with bright knots of star-formation in its spiral 
arms and with its shape possibly resembling NGC 1300. The ap- 
proximate position of the radio transient is shown by the circle 
(following Jakobsson et al. 2005). The field of view is 1 1 x 11 
arcsec^. North is up, east is left. 



3. Spectral energy distribution 

The spectral ene rgy distr ibution (SED) was analysed using 
HyperZ (Bolzonel la et al. 1 12000). The redshift of the host was 
fixed to the spectroscopic value of z = 0.41. We allowed the 
extinction value Ay to vary and fit the host extinction by us- 
ing 4 different redde ning laws: MW, LMC, SMC, and Calzetti 
dCal zetti et al. "2000^. To fit the stellar populations, we used two 
different stellar population model libraries: the GISSEL library 
of iBruzual & Chariot I (1 19931) and the stellar population models 
of iMaraston . (l2005h . The initial mass function of all models fol- 
lows the Salpeter law dSalpeter 111955b . 

The GISSEL library contains models with 8 different star- 
formation histories: an instantaneous burst, a constant star- 
forming system, and six exponentially decaying star-formation 
rate (SFR) models where SFR oc exp(-t/T), r being the charac- 



teristic lifetime. There are two sets of stellar population models 
in the GISSEL library: i) those with metallicities fixed to the so- 
lar value, and ii) those for which the metallicity evolve, where 
the evolution in metallicity of the stellar population is explicitly 
taken into accoun t so that there is a natura l link between age and 
mean metallicity (iBolzonella et al. Il2000l) . The ages of the tem- 
plates vary between zero and the age of the Universe (at z = 0.41 
in the case of th is GRB) covered with unevenly distributed steps 
(IBolzonella et al. 2000). 

From IMaraston i (i2005i) models, we used both single and 
composite stellar population models. We selected the models 
with metallicity Z - 0.04 (twice the solar value) to maintain 
agreement with the metallicity measured based on the spectro- 
scopic line ratio of rNIIU6584/rOIIU3 727, leading to log(0/H) 
-I- 12 = 9.0+0.1 d Levesque et al. 11201 Ol) . The single stellar pop- 
ulation (SSP) models that we used have a red horizontal branch 
morphology, meaning that the entire horizontal branch lifetime 
is spent on the red side of RR Lyrae strip (.Maraston ..2005.) . 

The host galaxy of GRB 020819 was also observed using 
IRAC on-board the Spitzer satellite. Svensson et al. I (1201 Ol) mea- 
sured the IRAC 3.6 fim and 5.8 jum brightnesses as 18.96+0.02 
mag (corre sponding to 95.1 + 1. 8 yuJy) and 19.27+0.22 mag, re- 
spectively. ICastro Ceron et al7l (l2008h also reports the IRAC 
3.6 //m band brightn ess as 97±2 fiJy, in agreement with 
'Svensson et al. I (|2010). We used these measurements with our 
GROND data to form a broad SED that covers the rest-frame UV 
to infrared wavelengths. In the following paragraphs, we present 
the SED fits of the host galaxy of GRB 020819, using GROND 
data combined with the IRAC brightnesses. The errors quoted 
in the following sections are at the 1 cr level and calculated by 
taking into account all the galaxy templates of a given stellar 
population library used in that fit. 

3.1. Fits to the GROND and IRAC data 

The best-fit models to the broad-range SED formed by GROND 
and IRAC data, are those with a MW-like extinction for both 
stellar library models (see Fig.|2]i. Using the models with evolv- 
ing metallicity from the GISSEL library, the best fit galaxy 
model has an age 30^j^'^ Myr and an extinction of Ay = 
2.21+0.39 mag using a MW reddening law and an exponen- 
tially decaying SFR with r=l Gyr (xl - 2.73). The stellar mass 
of the best-fit model i s -10''^'* solar masses, in agreement with 
ISavaglio et aTl d2009l) . 

Using Maraston SSP and CSP models with a fixed metallic- 
ity of Z = 0.04, the best-fit is obtained with a CSP model that has 
an exponentially decaying SFR with t=1 Gyr, age 6.0+0.5 Myr 
and Av=2.47;^Qj3 using a MW reddening law. The models with 
longer SFR e-folding times, r, provided similarly good fits with 
the same age and extinction values as best-fit. All other redden- 
ing laws provided a significantly (>1 cr) poorer fit. SSP models 
also fit similarly well but only when using a Calzetti reddening 
law with a best-fit model of Av=3.12 and age = 3 Myr 

We note that when using a MW reddening law, even the 
best fit does not provide a good fit to the H-band (see Fig. |2]i. 
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Fig. 2. The fits to the SED of the host of GRB 0208 19 using GROND and IRAC data. Note that we repeated the fits after increasing 
the IRAC 3.6/im band brightness error from 0.02 mag to 0. 1 mag. Error bars for both cases are shown in the figure. 



Similarly, the stellar population models extinguished with the 
Calzetti reddening law overestimate the /T-band flux. The rea- 
son for the imperfect fits to the //-band flux is possibly emis- 
sion lines such as Pa^ and He I from the ionised interstellar gas, 
which fall into the //-band at the redshift of this galaxy. 

4. Results and discussion 

The SED fits show that the optical and near-infrared emission 
from the galaxy is dominated by a very young stellar popula- 
tion. The young age of the dominant stellar population prevents 
us from distinguishing between models with different SFR e- 
folding times and hence does not allow us to determine the SFR 
history of the galaxy. The age of the dominant stellar population 
is ~6 Myr, if the metallicity is fixed to Z = 0.04, and slightly 
older, ~30 Myr, if the metallicity evolves in proportion to the 
age of the stellar population. The MW-like extinction provides 
better fits to the data than the Calzetti reddening. The best-fit ex- 
tinction value is in the range Ay ~ 1.8 - 2.6 mag depending on 
the galaxy template library used. 

In principle, there is a degeneracy between the dust extinc- 
tion in a galaxy and the age and metallicity of the mean stel- 
lar population of a galaxy, since all three redden the colour of 
the galaxy. However, it is also known that metal-poor galax- 
ies have less d ust content and lower global infr ared emission 
(ICalzetti et al. I l20IOt ICannon et al. 1 120051 l2006t IWalter et al. I 
l2007h . This is supported by our finding that the dust extinction 
and the age values obtained by fitting the solar/super-solar metal- 
licity models and the evolving metallicity models roughly agree 
with each other, indicating that either the metallicity is high or it 
does not have a significant eff'ect on the determination of the dust 
extinction for this galaxy, or both. In other words, the galaxy is 



highly extinguished according to the fits by synthetic stellar pop- 
ulation models, independent of the metallicity of the models. 

The age and dust extinction va lues we obtained agree with 
the spectroscopic measurements of iLevesque et al. I (|201C|) who 
found an age of 7.8+0.9 Myr using the Hyj emission line equiva- 
lent width and E (B - V) = 0.64 using the H„/H/3 line ratio. For 
a MW-like total-to-selective extinction ratio Ry = 3.1, E (B - V) 
- 0.64 corresponds to Ay - 1.98 mag, which is in agr e emen t 
with the Ay values from our SED fits. ILevesque et al. I (1201 Oh 
also measured the metallicity of the host galaxy based on the 
spectroscopic line ratio of [NII]/l6584/[OII]/i3727 and found a 
super-solar metaflicity of log(0/H) -H 12 = 9.0+0.1. 

The age (and extin ction) value obta i ned fr om either our SED 
fits or t he spectrum, bylLevesaue et al. I (l20IOh . do not agree with 
that of I Savaglio et al. I (l2009b . The reason for this disagreement 
is probably that they fit the SED of this host using only the B- 
, R-, and /T-band magnitudes available in the literature, where 
the B-band magnitude has a large uncertainty (see Sect. 2). 
Furthermore, their fits are restricted to the Ay range of 0.0 - 
2.0 mag and metallicities no more than the solar value, which do 
not agree with the best-fit model Ay or metallicity of this host. 
Hence, our analysis illustrates the importance of accurate pho- 
tometry, a broad SED coverage, and correctly selected parameter 
ranges. 

On the other hand, our bes t-fit stellar mass of ~ 10'°"^Mh 
agrees with ~1O'°-^"=''"'*M0 of SavagUo et al. I (l2009 l). indicat- 
ing that the stellar mass calculations based on the optical/near- 
infrared SED are much less affected by the choice of parameter 
ranges. A reason for the robustness of the mass estimation is 
that the SED covers the Balmer break and the rest-frame near- 
infrared emission of the galaxy. 
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The calculation of the dust extinction based on the after- 
glow emission limits of GRB 020819 in the R- and /T'-band by 
[Klose et al. (2003) show that if the afterglow is similar to that 
of GRB 011121 and at the same redshift (z=0.36), then an ex- 
tinction of Ay ~ 4 mag is required to dim the /?-band afterglow 
and Ay ~ 20 mag to dim the /T'-band afterglow below the lim- 
its. We performed a fit to the afterglow upper limits of the R- 
and /T'-band at z = 0.41 using a host galaxy extinction of Ay 
= 2.2 mag (following a MW-like extinction law) and found that 
the afterglow should have a spectral decay index/? >0.82. This 
value is very close to the canonical value of 0.8, showing that the 
afterglow of this dark GRB is probably quite normal i.e., not sig- 
nificantly different from other long GRB afterglows. Therefore, 
the "darkness" of this GRB is probably not related to the prop- 
erties of its afterglow but the dust in the host galaxy. The high 
Ay value is consistent with the darkness of the GRB afterglow, 
although the Ay that we deduced here is not only for the GRB 
region but the whole galaxy. 

There are several bursts with high extinction in the circum- 
burst region or along the line of sight, derived using the af- 
terglow emission (e.g. Diorgovs ki et al. 2 001; Rol et al. 2007; 
Jaunsen et al. 2008; Ta nvir et al. 2008; K ann et al. I2006ll2009t 
Hashimoto et al. „2010i) . However, GRB 020819 is the first case 



with known redshift and a well-constrained high extinction value 
for the entire galaxy. A similar case is that of GRB 0904 17B 
(z=0.345) for which Holland et al. (2010) measure Ay -3.5 
mag based on the UV and optical SED of this dark GRB 
host. The high host-galaxy extinctions of these low -redshift dark 
GRBs also show that while the dust extinction is one of the main 
reasons for the darkness of low-redshift bursts, at high-redshift 
(z>l) the dust extinction required to dim the afterglow is not that 
high. 

The case of GRB 020819 is an exception in the context of 
both the host galaxies of long GRBs with bright optical after- 
glows, and the dark GRB hosts. Host galaxy studies of dark 
GRBs indicate that the majority of the dark GRB hosts are sim- 
ilar to normal long GRB hosts , which do not have very high ex- 
tinction or very high redshifts (iPerlev et al. 1 2009). However, the 
host extinction of GRB 020819 is much higher than the average 
extinction. Ay < 0.5, of a sample of hosts galaxies of long GRBs 
with detected optical afterglo ws, obtained either using spec- 
trosco py or SED fitting (e.g., iHan et al.ni2010t ISavaglio et al. I 
12009 ^. Furthermore, the spectroscopy of the host indicat es very 
high star-formation-rate (SFR) of 23.6 Mo/yr (Leves aueet al. I 
1201 Ol) corresponding to a specific SFR (SFR / stellar mass) of 
0.94 per Gyr ' using our best-fit stellar mass IO'^'^Mq. This spe- 
cific SFR is higher than the median specific S FR ~0.8 Gyr"' of 
long GRB host galaxies (Savaglio et al. ''2009'). These properties 
mark the host galaxy of GRB 020819 as a galaxy with dusty, 
intense star-formation, distinguished from the currently known 
GRB host population. 

The host galaxy of another dark GRB 051022 (z=0.807) 
shows similarly exceptional properties in terms of mass (~ 
lO'^'^Mo), S FR (~ 36 M»/yr) and metallicity ( l og(0/ H) + 
12 = 8.77) (ISavaglio et al. I (20091: iGraham et al.1 l2009h . We 
note however that there is a large uncertainty in the metal- 
licity calculation for this host depending on whether the up- 
per branch solution calibrated with metal-rich galaxies or a 
lower branch solution calibrated with metal-poor galaxies i s 
used (see e.g.. lKobulnickv et al. ll999l:lKewlev & Dopital2002h . 
The extinction of the host galaxy of GRB 051022 (Ay=1.55, 
iLeve sQue et al. 2010) is also relatively high compared to the av- 
erage long-GRB host population. Similarly, there are a few dark 
GRBs that have both high extinction and near-solar metallicities 



derived using aftergl o w observations (e.g.. Watson et al. II200l 
'Kruhler et all l2008t lEliasdottir et al. 1 120091: iProchaska et al 
2009), indicating that hosts such as that of GRB 020819 may 
be more common than previously anticipated. 

These findings and the clear example of the host galaxy of 
GRB 020819, show that the hosts of GRBs may be far more 
diverse, in terms of extinction, mass and metallicity, than previ- 
ously found. 
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